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Abstract 

We provide a unified theory of luminescence spectra of coupled light-matter systems realized 
with semiconductor heterostructures in microcavities, encompassing: i) the spontaneous emission 
case, where the system decays from a prepared (typically pure) initial state, and ii) luminescence 
in the presence of a continuous, incoherent pump. While the former case has been amply discussed 
in the literature (albeit mainly for the case of resonance), no consideration has been given to the 
influence of the incoherent pump. We show how, by provoking a self-consistent quantum state, the 
pump considerably alters the emission spectra, even at vanishing intensities. The main outcome of 
our analysis is to unambiguously identify strong-coupling in situations where it appears in disguise 
or only seems to appear. Here, we consider bosonic matter fields, in which case fully analytical 
solutions can be obtained. This describes the case of quantum wells or large quantum dots, or the 
limit of low excitation where the average populations remain much smaller than one. 

PACS numbers: 42.50.Ct, 78.67.Hc, 42.55. Sa, 32.70.Jz 
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I. INTRODUCTION 



The dynamics of an optical emitter changes drastically when it is placed in a cavity. 
The cavity alters the density of states of optical modes, and therefore increases or inhibits 
interactions with the emitter. The effect was first put to use by Purcell in nuclear magnetic 
resonance for the practical purpose of thermalizing spins at radio frequencies, by bringing 
down their relaxation time from ^ lO^'^s to a few minutes Kleppner applied the same 
idea in the opposite way, to increase the relaxation time of an excited atom, i.e., to inhibit 
its spontaneous emission (SE) [2f|. The emitter, that in the case of Purcell was sought to 
be resonant with the cavity mode to increase the photon density of states with respect 
to the vacuum, was in the case of Kleppner put out of resonance, namely in a photonic 
gap, where the photon density of states is smaller than in vacuum. This tuning of the 
relaxation time of an emitter placed in a cavity, now known as the Purcell effect, has many 
potential technological applications, one of the most compelling being the decrease of the 
lasing threshold. The effect, that had first been actively looked for with atoms in cavities j^, 
was therefore also intensively (and more recently) pursued in the solid state, more prone 
for massive technological implementations. Semiconductor heterostructures are the state of 
the art arena for this purpose. They allow to engineer, with an ever rising control, the solid 
state counterpart of the atomic system to match or isolate their excitation spectra and thus 
control their behaviour. Typical examples are quantum dots (QDs) placed in cavities made 
in micropillars, microdisks or photonic crystals, where Purcell inhibition has been distinctly 
demonstrated 4, l5|. In this paper, we shall be concerned chiefly with the observed spectra 
of emission of such an emitter placed in a cavity. In the regime where the effect of the cavity 
is to lengthen or shorten the lifetime of the excitation, the consequence in the optical spectra 
is to narrow or broaden the line, respectively. 

In the description of the Purcell effect, the possible reabsorption of the photon by the 
emitter is so weak that it can be neglected. It is responsible for the energy shift known as 
the Lamb shift, that, in quantum electrodynamics, is interpreted as the perturbative influ- 
ence of virtual photons emitted and re-absorbed by the emitter. In the case of inhibition of 
the spontaneous emission, this shift is indeed orders of magnitude smaller than the radia- 
tive broadening. In the case where emission is enhanced, and the linewidth narrowed, the 
probability of reabsorption of a photon by the emitter becomes closer to that of escaping 
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the cavity, until the perturbative — so-called weak-coupling (WC) — regime breaks down and 
instead strong coupling (SC) takes place. In this case, photons emitted are then reflected by 
the mirrors and there is a higher probability for their reabsorption by the atom than for their 
leaking out of the cavity. A whole sequence of absorptions and emissions can therefore take 
place, known as Rabi oscillations. This regime is of greater interest, as it gives rise to new 
quantum states of the light-matter coupled system, usually referred to as dressed states in 
atomic physics and as polaritons in solid-state physics. Experimentally, SC is more difficult 
to reach, as it requires a fine control of the quantum coupling between the bare modes and 
in particular to reduce as much as possible all the sources of dissipation. Theoretically, it is 
better dealt with by flrst getting rid of the dissipation, and starting with the strong- coupling 
Hamiltonian {h is taken as 1 along the paper): 

+ LUbb^b + g{a^b + ab^) (1) 

where a and b are the cavity photon and material excitation fleld operators, respectively, 
with bare mode energies Ua and Uh, coupled linearly with strength g. The photon operator is 
a Bose annihilation operator, satisfying the usual commutation rule [a, a^] = 1. Depending 
on the model for the material excitation, b is described by, typically, another harmonic 
oscillator (linear model Q|) or a two- level system (Jaynes-Cummings model j?!). Those are 
the most fundamental CclSGS cLS they describe material flelds with Bose and Fermi statistics, 
respectively. Possible extensions are a collection of harmonic oscillators jsl or of two-level 



systems [9|, a three-level system 10|], etc. This paper contains the flrst part of our work 
in which we address exclusively the case where b also follows Bose statistics. This is an 
important case for two reasons. The flrst one is that in many relevant cases, the matter-fleld 
is indeed bosonic, such as the case of quantum wells, or large quantum dots, at low density of 
excitations. The second reason is that this case provides the limit for vanishing excitations of 
all the other cases, and is fully solvable analytically. As such, this case provides the backbone 
for the qualitative understanding of SC. We leave for a further work the case of fermions at 
larger pumping, more relevant when dealing with small QDs that conflne excitations, and 
more prone to involve genuine quantum mechanics as one quantum of excitation can alter 
the system's response. The drawback is that numerical computation is required in this case, 
and the discussion is therefore of a somehow less fundamental character. 

For the rest of the text, it is therefore understood that b is also a Bose operator. For 
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convenience, we shall refer to it as the exciton operator, after the name of a bound electron- 
hole pair in a semiconductor. Likewise, we shall prefer such terminology as a Quantum Dot 
(QD) rather than an atom, or polaritons, rather than dressed states, etc. For most purpose, 
this is semantics only and the results apply in a wide range of systems. 

Equation ([T]) can be straightforwardly diagonalized for bosonic modes, giving H = 
uJvP^P + ^^hQ^Q, where 

-u = ^— ±^^^ + -, (2) 
with new Bose operators p = {cos9)a + (sin6')6 and q = — (sin^)a + {cos6)b, determined by 



the mixing angle, 6 = arctan ((y + yS'^ + 'x)/^')' ^^'^ detuning: 

A = Ua - Ulb ■ (3) 

These new modes are the polaritons (or dressed states) with quantum states |U) = |vac) = 
cos^ |1, 0) + sin^ |0, 1) and |L) = g"!" |vac) = — sin6' |1, 0) + cos^ |0, 1), where |vac) is the 
vacuum, |1, 0) is the Fock state of one photon and |0, 1) the Fock state of one exciton. 

The energies defined by eqn ([2]) are displayed in Fig. [1] with dashed lines, on top of that of 
the bare modes, with thick lines, as detuning is varied by changing the energy of the emitter 
and keeping that of the cavity constant. The anticrossing always keeps the upper mode U 
higher in energy than the lower L one, strongly admixing the light and matter character of 
both particles. At resonance, this mixing is maximal: 

IE) = ^(|1,0)±|0,1)) . (4) 

If the system is initially prepared as a bare state — which is the natural picture when reaching 
the SC from the excited state of an emitter — the dynamics is that of an oscillatory transfer 
of energy between light and matter. In an empty cavity, the time evolution of a state which 
was an exciton at t = 0, is given by: 

= exp ( - iHt) |0, 1) , (5a) 
= sin^cos^(e-*'^u* -e-^'^Lt) |i^o) 

+ (sin^ ^e-^'^ut + cos2 ^g-^^^L*) |o^ i) ^ (5b) 
which results in oscillations between the bare modes at the so-called Rahi frequency: 
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FIG. 1: (Color online) Solid black: Bare energies of the cavity photon (horizontal line) and of 
the exciton (tilted) as a function of detuning A. Dashed black: Eigenenergies of the system 
Hamiltonian, without dissipation nor pumping [eqn ([2|)]. The exciton-like state at large negative 
A has become a photon-like state at large positive A, and vice-versa. Around A = 0, both modes 
are an admixture of exciton and photon. Dotted blue: Correction of the eigenenergies when pump 
and decay are taken into account for parameters of point (c) in Fig. [9l Solid red: Actual position 
of the observed peaks in the PL spectra for the same parameters. For these parameters, the three 
descriptions of SC give the same qualitative results. 

The emission spectrum of the system requires dissipation, as it is an obvious practical 
requisite that the excitation should eventually leak out of the system to be detected from 
the outside [the energy spectrum of the system is given by eqn ([2])]. Dissipation is intrinsic to 
the bare modes: both the cavity photon and the exciton have a finite lifetime. In presence of 
dissipation, the system is upgraded from a Hamiltonian description, eqn ([1]), to a Liouvillian 
description, with a quantum dissipative master equation for the density matrix p defined in 



the tensor product of the light and matter Hilbert spaces Tia and Tib 
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dtp = C^p . (7) 

The Liouvillian C-y still contains the Hamiltonian dynamics of SC, but also takes into account 
the decays of both the cavity and the emitter, with rates 7a and 7^, respectively: 

C^p = i[p, H]+Y^ jC^cpc^ - c^cp - pc^c) , (8) 



in which we have considered temperature equal to zero. This equation has been extensively 



studied 



12l |. although, to the best of our knowledge, not in its most general form. The 



typical restrictions have been to consider the case of resonance, uja = uJb, with only one 
particular initial condition, namely, the excited state of the emitter in an empty cavity, and 
to detect the emission of the emitter itself. All together, they describe the spontaneous 
emission of an emitter placed into a cavity with which it enters into SC. This has been the 
topical case for decades as this was the case of experimental interest with atoms in cavities. 

With the advent of SC in other systems, other configurations start to be of interest. 
With a QD in a microcavity, the detuning A between the modes, eqn ([3]), is a crucial 
experimental parameter, as it can be easily tuned and to a great extent, for instance by 
applying a magnetic field or changing the temperature. Also in this case, the detection is 
in the optical mode of the cavity, rather than the direct emission of the exciton emission, 
because the latter is awkward for various technical reasons of a more or less fundamental 
character (an example of a fundamental complication is that the emission is enhanced in the 
cavity mode and suppressed otherwise, and the exciton lifetime is typically much longer, so 
the exciton emission is much weaker; an example of a petty technical complication is that the 
exciton detection should be made at an angle and, practically, a lot of samples are grown 
on the same substrate. Both the substrates and other samples hinder the lateral access 
to one given sample, whereas all are equally accessible from above). In our system where 
both modes are bosonic, symmetry allows us to focus on the cavity emission without loss of 
generality, as we can obtain the leaky excitonic emission by simply exchanging indexes a, b 
(the spectrum could also have photon-exciton crossed terms that could be computed in a 
similar way). When we shall turn to the case of a fermionic matter-field, where the exciton 
emission will become distinctly different and for that reason, important, we shall address 
exciton emission separately 13 1. 

Regarding the initial condition, more general quantum states can now be realized, at 
least in principle, by coherent control, pulse shaping or similar techniques. Additionally and 
more importantly, the type of excitation of a cavity-emitter system in a semiconductor is of 
a different character: the excitation is typically injected by either a continuous wave (cw) 
laser far above resonance, creating electron-hole pairs that relax incoherently to excite the 
QD in a continuous flow of excitations, or by electrical pumping as in lasers. This pumping, 
that is of an incoherent nature typical of semiconductor physics, brings many fundamental 
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changes into the problem that go beyond the mere generahsation of eqn ([8]). Among the most 
obvious ones, let us already mention that pure states of the like of eqn (jlj) do not correspond 
to the experimental reality. Instead, the system is maintained in a partially mixed state with 
probabilities p{n) to realize the nth excited state. In all steady state is imposed by 

the interplay of pumping and decay. The Rabi oscillations of the populations — that is, the 
coherent exchange of energy between the modes — are always washed out, regardless of the 
photon-like, exciton-like or polariton-like (eigenstate) character of the density matrix. 

In this work, we address both the emission spectra obtained in a configuration of spon- 
taneous emission (SE) — where an initial state is prepared and left to decay, as ruled by 
eqns ([THH]) — under its most general setting, and the case of luminescence emission under 
the action of a continuous and incoherent pumping that establishes a steady state (SS). We 
bring all results under a common and unified formalism and show how none of the cases 
fully encompasses the other. The rest of the paper is organised as follows. In Section [Tll 
we present the complete model and we derive and discuss its equation of motion, focusing 
on the single-time dynamics. In Section IIIIl we obtain fully analytically the main results 
in both of the cases explicated above, this time focusing more on the two-time dynamics, 
which Fourier transform gives the luminescence spectra. In Section [TVl we discuss the math- 
ematical results derived in the two previous sections, accentuating the physical picture and 
relying on particular cases for illustration. We consider, in this Section, specifically the case 
of resonance, where all the concepts manifest more clearly. Finally, in Section |Vl we give a 
summary of the main results and provide an index of all the important formulas and key 
figures of this text. We conclude with a short overview of the continuation of this work that 
replaces the bosonic emitter with a fermionic one. 

II. MODEL 

Our model for the coupling of a QD with a single cavity mode in the presence of incoherent 
and continuous pumping is sketched in Fig. [2l One QD is strongly coupled to the cavity mode 
with interaction strength g, and is continuously, and incoherently, excited by an electronic 
pumping Pb, which, in the microscopic picture, is linked to the rate at which electron- 
hole pairs relax into the dot. This rate is related in some way to the pumping exerted 
by the experimentalist. Because of the incoherent nature of this pumping, together with 
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FIG. 2: (Color online) Schema of our system for the SS case: self-assembled QDs in a semiconductor 
microcavity. The QD sketched on the right is in SC with the cavity mode with coupling strength g, 
while the one of the left is in WC. The electron- hole pairs created by the incoherent pumping of 
the structure provide an effective electronic pumping, Pf,, of the SC dot, while the pumping of the 
WC dot results in an effective cavity pumping Pa, through rapid conversion of the excitons into 
cavity photons. 

that of the damping, the off-diagonal elements of the QD reduced density matrix (that hold 
the coherence) are washed out in the SS. We also consider another type of pumping. Pa, 
that offers a counterpart for the cavity by injecting photons incoherently at this rate. The 
major factor to account for such a term is the presence of many other QDs, that have been 
grown along with the one of interest. Those only interact weakly with the cavity. In most 
experimental situations so far, it is indeed difficult to find one dot with a sufficient coupling 
to enter the non-perturbative regime. When this is the case, all the other dots that remain 
in WC become spectators of the SC physics between the interesting dot and the cavity, 
and their presence is noticed by weak emission lines in the luminescence spectrum and an 
increased cavity emission. They are also excited by the electronic pumping that is imposed 
by the experimentalist, but instead of undergoing SC, they relax their energy into the cavity 
by Purcell enhancement or inhibition, depending on their proximity with the cavity mode. 
This, in turn, results in an effective pumping of the cavity [l^ . 

To model these two continuous and incoherent pumping, we use Lindblad terms Cp 
that substitute the annihilation operator c (standing for a or b) by the creation one c^, 
and vice-versa, a procedure which is known to describe pumping terms in a quantum rate 
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equation [15|]. The full Liouvillian now reads (again at zero temperature): 

dtp =Cp = {C, + Cp)p (9a) 
=iuJa[p, a^a] + i^bip^ ^^^] + ig[p-i o-^b + ab^] (9b) 
+ ^{2cpc'^ - c^cp - pc^c) (9c) 

p 

+ ^ y(2cVc - ccV - pcc^) ■ (9d) 

c=a,6 



The microscopic derivation of line (I9dl) follows from the usual Born-Markov approxima- 



16[ . Since 



tion ll|. The case of electronic pumping, for instance, is similar to the process of laser gain: 
the medium requires an inversion of electron-hole population, something that cannot be 
achieved by means of a simple harmonic oscillator heat bath. The actual process of gaining 
an exciton in the QD involves the annihilation an electron-hole pair in an external reservoir 
out of equilibrium (representing either electrical injection or the capture of excitons optically 
created at frequencies larger than the ones of our system) and the emission of a phonon, that 
carries the excess of energy, to another one (which can be in thermal equilibrium). A simple 
effective description of this non-equilibrium process can be made by an inverted harmonic 
oscillator with levels Ep = —up{p+ 1/2) maintained at a negative temperature 
the raising operator for the energy decreases the number of quanta of this oscillator, the 
role of creation and destruction operators is indeed reversed with respect to the usual case 
of damping. 

For the sake of generality, we introduce effective broadenings, that reduce to the decay 
rates in the SE case but get renormalized by the pumping rate in the SS case: 

Ta.fe = 7a,fe , (SE case) (10a) 

^a,b = la,b - Pa,b ■ (SS CaSc) (10b) 

We shall also use thoroughly the combinations: 

7± = — ^ — and r± = ^ . (11) 

The renormalization of the linewidth in the presence of the pumping term is a bosonic 
effect. In the case of a single harmonic oscillator driven by pump and decay, the emission 
spectrum is a Lorentzian lineshape with Full- Width at Half Maximum (FWHM) given by 
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7 — -P = P/ {n). The linewidth narrows with the number of particle, in a way reminiscent 



of the Schallow-Townes effect 



17|. 



From the relations (O) = Tr(Op) and dt{0) = Ti{Odtp) = Tr(0£p), we can obtain from 
eqn (jH]) the single-time mean values of interest for this problem, by solving the equation of 
motion of the coupled system: 



(12) 
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and providing the initial conditions: 

^l^^ai^), n|J = rai,(0), and n^^ = ra„b(0) . 



(13) 



The solutions are heavy, but not devoid of interest. For completeness, and since this is a 
natural and potentially useful result within the scope of this paper, we provide them in 
Appendix \^ 

On the other hand, the SS case corresponds to setting the time derivative on the left 
hand side of eqn f|T2l) to zero, and solving the resulting set of linear equations. This yields: 

^^_g^T^{P, + Pb)+PaTb{Tl+{^f) 



n 



Ag^Vl + VJ^biXl + (1)2) 



SS 



„ss 

^ab 



^2r+(p, + n) + Ar,(r^ + (f)2) 

Ag^Tl + Tj:b{Tl + (f )2) 
f(7an-7fePa)(^r+-f) 



(14a) 
(14b) 
(14c) 



4^72r^ + r„r,(r^ 

Both photonic and excitonic reduced density matrices are diagonal. They correspond to 



thermal distributions of particles with the above mean numbers 

„ ,,,SS^r^ 



m; p,m 



Pm,q 



m 

' ^ ^ Pn,m; n,q 



in': 



n,p 



1 + nf )"+! 



in' 



SS\m 



m,q 



1 + np)"^+i ■ 



(15a) 
(15b) 



Behind their forbidding appearance, eqns enjoy a transparent physical meaning, that 
they inherit from the semi-classical — and therefore intuitive — picture of rate equations. For 
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instance, when the couphng strength between the two modes, g, vanishes, the solutions 
na,b = Pa,b/^a,b are those of a source and sink problem for bosons, i.e., of the kind dtn = 
—'yn + P{n + 1), featuring the famous Bose stimulation effect, whereby the probability of 
relaxation towards the final state is increased by its population. 



III. CORRELATION FUNCTIONS AND SPECTRA 

We now turn to the main goal of this paper, namely, the luminescence spectrum of the 
system s(u;). Physically, it is the intensity (or mean number) of photons in the system with 
frequency u, i.e., by definition: 

s{uj) = {a^{uj)a{uj)) . (16) 

This is proportional to the intensity of photons emitted by the cavity at this frequency (the 
direct exciton emission from its recombination is described in a similar way by {l)^ {uj)h{uj))) . 
It will be more convenient, throughout, to deal with normalized spectra: 

S{lo) = s{u) I {a)a){t)dt, (17) 

so that eqn f|T7|) is now the density of probability that a photon emitted by the system 
has frequency uo. The Fourier transform of a{ijj) relates the emission spectrum to a two- 
time correlator through S{uj) = ^ JJ^ {a'^ {ti)a{t2))e'^^^*^~^^^dtidt2/ f {a^a){t)dt, that after a 
change of variables, can be expressed in terms of the so-called first-order time autocorrelator: 

, r) = {a\t)a{t + r)) (18) 

for positive time delay r = t2 — > 0. All put together, this yields the usual Fourier-pair 
relationship between the power spectrum and the autocorrelation function: 

= ^^^^TTTTTT:^^ / / G^^\t,T)e'-^drdt. (19) 
TT jg {a^a)[t)dt Jo Jq 

Equation (fT9l) holds as such in the SE case. In the SS case, care must be taken with 
cancellation of infinities brought by the ever-increasing time t. A technical but straight- 
forward procedure (cf. Appendix [B]) leads to the expression that explicitly gets rid of the 
divergences — famously known as the Wiener- Khint chine theorem |l^ — that reads: 

S^^{u) = 5^ lim 3f? / G(^)(t,r)e'""rfr. (20) 
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From now on, we shall refer with "SE" and "SS" to the expressions that apply specifically 
to the spontaneous emission and to the steady state, respectively, leaving free of index those 
that are of general validity. In some cases, as for instance in eqn fllOp . no index is required 
if it is understood that Pa/b are defined and equal to zero in the SE case. For that reason, 
we shall leave F free of the SE/SS redundant index. 

To obtain the spectra of a system whose dynamics is dictated by eqn (j9]), we therefore need 
to compute two-time dynamics. This can be done thanks to the quantum regression theorem, 
according to which, a set of operators Aj^^} that satisfy Tr(A{Q,}£0) = X]{/3} ^{a/3}Tr(^{/3}0) 
for all O G Tia ® 'Hb for some M{q^}, yields the equations of motion for the two-time 
correlators as: 

dr{0{t)A{^^{t + r)) = Y, M{,^}(0(t)A|,|(t + r)) . (21) 

m 

In the fully bosonic case, operators a"6™ with {a} = (n, m) G N constitute such a set 
with M defined by 

p 

Mnm = —iinuja + mujb) — n— — m — , (22) 

nm 2 2 

M nm = M mn = —iom , (23) 
n+l,m — 1 m— l,n+l 

and zero everywhere else. For the computation of the optical spectrum, it is enough to 
consider = a and O = a''' in eqn ( 1211) . We obtain the equation 



for the correlators 



a,v(t,t + r) = Miv(t,t + r) (24) 



where 



v(M + r)H^'''<"''""'"l (25) 

(at(t)6(t + r)) 

Mio Mio \ l-iuja -ig ^ 

Moi Moi / \ -ig -iuJb 

10 01/ \ 

The formal solution follows straightforwardly from v(t, t+r) = e'^^'^v(t, t). Made explicit, 
it reads (at positive r): 

exp[-(F+ - 'if )r] 



{a\t)a{t + T)) 



2R 

X |e^^^[(i? + tT_ - A/2)na{t) -gnab{t)] 

+ e-^^^[(i? - zF_ + A/2)na{t) +gnab{t)]} (27) 
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in terms of the complex Rabi frequency: 

R = 



.A 



(28) 



that arises as a direct extension of the dissipationless case, eqn For our discussion, 
it is convenient to decompose R into its real and imaginary parts, R = R^ + iR\. Out of 
resonance, the Rabi frequency is a complex number with both nonzero real and imaginary 
parts. On the other hand, at resonance, it is either pure imaginary (in the WC regime), 
either pure real (in the SC one). For this latter case it is worth defining a new quantity: 



R, = R{A = 0) = Jg^ - Tl . 



(29) 



In Fig. [3|, the real and imaginary parts of R are plotted as a function of T^/g for various 
detunings. In the limit of high detuning, |A| ^ g, regardless of WC or SC, the real part 
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FIG. 3: (Color online) Complex Rabi R/g, separated in its real (a) and imaginary (b) parts, as 
a function of the decoherence parameter T^/g for various detunings {A/g from —1.6, up, to 0, 
bottom, by steps of 0.4). Solid black lines correspond to resonance. 

becomes independent of the dissipation (decay and pumping), Rj. ^ | A|/2, and the imaginary 
part becomes -Ri ~ ^Ff-- Therefore, as already pointed out, the limit of bare modes at 
energies Ua = and Ub = —A broadened with the bare parameters Ta/b (FWHM), is always 
recovered at large detunings. 

Using the result of eqn (1271) into the definition of eqn (1191) . we obtain the formal structure 
of the emission spectrum: 



S{uj) = — {C^ + C^) - ^{C}{C^ - - ^{C]{A^ - 
2tx 



(30) 
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with and A{uj) some Lorentzian and dispersive functions whose features (position and 
broadening) are entirely specified by the complex Rabi frequency [eqn fl28|) ]. r+ [eqn ffTTl) ] 
and the detuning A: 



^'''(^) =7^^^7T^^ ; A^o^^. ' (31a) 



,1,2/ N ^ - (-T =F ^r) 



We also introduced C, a complex coefficient given by 



C = R ^ (32) 



that we define in terms of still another parameter, D: 



{a%){t) dt 

D = ^ . (33) 

{a)a){t)dt 



Written in this form, eqns (jnOHSSD assume a transparent physical meaning with a clear 
origin for each term. The spectrum consists of two peaks (that we label 1 and 2), as is 
well known qualitatively for the SC regime. These are composed of a Lorentzian C and 
a dispersive A part. The Lorentzian is the fundamental lineshape for a system with a 
lifetime, and in the expression above, it inherits most of how the dissipation gets distributed 
in the coupled system, including the so-called suhnatural linewidth averaging that sees the 



broadening at resonance below the cavity mode width [12|. The dispersive part originates 
from the coupling as in the Lorentz (driven) oscillator. In our system, it stems from the 
driving of one mode by the other, because of the coupling. This decomposition of each peak 
in such terms is therefore entirely clear and expected. More quantitatively, the first peak, 
(e.g.,) is centred at — ^ — Rr and broadened by r+ + Ri. As Rj- > 0, this peak corresponds 
to the lower branch "L". 

So far, all the results hold for both cases of SE and SS. This shows that the qualitative 
depiction of SC is robust. This made it possible to pursue it in a given experimental system 
with the parameters of the theoretical models fit for another. This has indeed been the 
situation with semiconductor results explained in terms of the formalism built for atomic 
systems. 



14 



To be complete, the solution now only requires the boundary conditions that are given 
by the quantum state of the system. They will affect the parameter eqn fl33l) . that is 
therefore the bridging parameter between the two cases. In the next two sections, we address 
the two cases and their specificities. 



A. Case of Spontaneous Emission 



In the case of Spontaneous Emission 12|, |20|, where the system decays from an initial 



state, the boundary conditions are supplied for r = by the initial values v(t, t), i.e., the cav- 
ity population, naif) = {a^a)(t) and the coherence element Uabit) = {a^b)(t). In turn, those 
are completely defined by the initial conditions, eqns (fT3|) . Although the analytical expres- 
sion for these mean values as a function of time are cumbersome — see Appendix \M — the D 
coefficient, eqn ( |33l) . that determines quantitatively the lineshape, assumes a (relatively) 
simpler expression: 

To prepare the analogy with the SS case in the next section, we also write the particular 
case when rfi,. = 0: 



f(7a<-7.n°)(z7+-f) 
^?27+(nO + nO) + n07,(72 +(fP) 



7-)SE _ 2^ '""ft i>J-~a/\" 2 ^ /qc\ 



This is an important case as it is realized whenever the initial population of one of the 
modes is zero. Note that in this case, D^^, and therefore also the normalised spectra, do 
not depend on the magnitude of the populations. 



B. Case of continuous, incoherent pumping 

In the case where the system is excited by a continuous, incoherent pumping, a steady 
state is reached and the boundary conditions are given by the stationary limit, as time 
tends to infinity, of the dynamical equation (whose solution is unique). The D parameter, 
eqn fl33|) . is defined in this 



nSS _ ''^ab _ 2V^"-'0 /o-^ an"--- + 2^ /nf:\ 



ss 



f(7an-7fePa)(^rH.-f; 



g^T^iPa + Pb) + PaTbin + (f )2) ■ 



15 



There is a clear analogy between eqn (l36l) — that corresponds to the SS — and eqn (l35i) — 
that corresponds to SE when ?7,°^ = 0. It is made more meaningful by defining the ratio 

a = n^/n'^ (SE case with ?7,°^ = 0) , (37a) 
a = Pa/ Ph (SS case) , (37b) 

in which case eqns (135|) and (136!) assume the same expression, keeping in mind the definition 
of eqns ffTOj) . Table [T] displays this common expression for D in terms of a. The limiting 
cases when a ^ or oo are also given. They correspond to only photons or excitons as the 
initial state for the SE, or to the presence of only one kind of incoherent pumping for the 
SS case. 

The analogy and differences between D^^ and D^^ reflect in the spectra S^^ and S"^^. 
For the same a, they become identical when the pumping rates are negligible as compared 
to the decays, P^^fe <^ '-)a,b- In this case, where T±^a,h ~ 7±,a,b5 the SS system indeed behaves 
like that of the SE of particles that decay independently and that are, at each emission, 
either a photon or an exciton, with probabilities in the ratio a. 

However, in the most general case, D^^ depends on more parameters — namely, a and 
r±,a,b — than _D^^, that depends only on a and 7±,a,fe, cf. eqn ffTOl) . Moreover, the pumping 
rates Pa^ affect S*^^ not only through a and D^^, but also in the position and broadening of 
the peaks (given by F-i- and R). Therefore, the SS is a more general case, from which the SE 
with = can be obtained, but not the other way around. On the other hand, as seen in 
Table [H the SS case cannot recover the SE case when ^ (this could be overcome if cross 
Lindblad pumping terms were introduced in eqn ^ with parameters Pab corresponding to 



the cross initial mean value n°{,, as is done in, e.g., Ref. 2l|, but this describes another 
system) . 

In any case, an important fact for the semiconductor community is that a SS with non- 



vanishing pumping rates is out of reach o 
the case studied in the literature so far 12 



;he SE of any initial state, which has been 



20|, and that even in this limiting case, the 



effective quantum state obtained in the SS should still be resolved self-consistently, rather 
than assuming for a the particular case or oo. 
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< Pa 


U 




< a < cx) 
oo 


-f (^r+ - f )7fe 
52rH_ + r,(r^ + (A)2)) 

f(ir+ - 4)(7a-7b") 


ff2r+(i + a) + arfe(r^ + (^)2)) 

(^r+ - t)7a 

2.9r+ 



TABLE I: Expression of D, eqn ([33]) . as a function of a, eqns ([37|) . in the SE (with T-i-^a^fe — > 7±,a,6 
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and = 0) and SS cases. D embodies in the luminescence spectrum the influence of the quantum 



state of the system. The latter is specified by the initial condition in SE, or the pumping/decay 
interplay in the SS. 

C. Discussion 



With this exposition of the analytical expressions of the luminescence spectra, and the 
discussion of their similarity and distinctions that we have just given, the coverage of the 
problem is complete and comprehensive. In order to give a more physical picture of these 
abstract results, we shall in the rest of this paper illustrate the implications that this bears 
in practical terms. We consider the case of resonance for this purpose, for reasons detailed in 
the next section. This will also allow us to provide self-contained expressions for the spectra. 



that can be used more conveniently to p 



fitting of the experimental data 
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ot the expression or use the model in a nonlinear 



24| . The system lends itself naturally to a global 



fitting, i.e., constraining all fitting parameters over the various detuning cases, where only 
the detuning is allowed to vary, while optimising them globally [l^. There is no difficulty in 
extending all the discussions that follow to arbitrary detunings. For instance. Fig. H] shows 
the SS spectra and their mathematical decompositions into Lorentzian and dispersive parts, 
as detuning is varied. Figs, (b) and (c) are obtained using eqn (IHOHSSl) . and in this particular 
case, the expression (l36l) for D. Fig. (a), the case at resonance, can be obtained in the same 
way, but in next Section we shall bring all these results together in a condensed expression 
[eqn 

More importantly, as we shall soon appreciate, the resonant case is the pillar of the SC 
physics. The main output of the out-of-resonance case is to help identify or to characterise 
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FIG. 4: (Color online) Strong-coupling SS spectra (blue solid line) and their decomposition into 
Lorentzian (dotted purple) and dispersive parts (dashed green) for various detunings (A/(7 = 0, 1, 3) 
with parameters of point (c) of Fig. (9) 7a = 3.8(7, Ih = O.lg, Pa = 0.5(7, Pb = 0.1(7. The vertical 
black lines mark the positions of the bare modes, showing the "level repulsion" of SC. 

the resonance, for instance by localising it in an anticrossing or by providing useful additional 
constrains with only one more free parameter in a global fitting. Even a slight detuning 
brings features of WC into the SC system and ultimately, when A ^ g, the complex Rabi 
frequency converges into the same expression for both regimes (as showed in Fig. [3]). This 
is why we now consider the SC problem in its purest form: when the coupling between the 
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modes is optimum. 



IV. STRONG & WEAK COUPLING AT RESONANCE 

Strong-coupling is most marked at resonance, and this is where its signature is experi- 
mentally ascertained, in the form of an anticrossing. Fundamentally, there is another reason 
why resonance stands out as predominant: this is where a criterion for SC can be defined 
unambiguously in presence of dissipation: 

WC and SC are formally defined as the regime where the complex Rabi frequency at 
resonance, eqn is pure imaginary (WC) or real (SC). 

This definition, that takes into account dissipation and pumping, generalises the classi- 
fication found in the literature, and is one of the main results of this work. The reason for 
this definition is mainly to be found in the behaviour of the time autocorrelator, eqn (l27j) . 
that is respectively damped or oscillatory as a result. The exponential damping is the usual 
manifestation of dissipation, that decays the correlations in the field, even when a steady 
state is maintained. On the other hand, in the same situation of steady averages (no dy- 
namics) but now in SC, the oscillations with r are the mark of a coherent exchange between 
the bare fields (the photon field and exciton field). 

In the literature, one sometimes encounters the confusion that SC is linked to a periodic 
transfer of energy or of population between the photon and exciton field, or that it follows 
from a chain of emissions and absorptions. This is an incorrect general association as one can 
explicit cases with apparent oscillations of populations that correspond to weak coupling, 
or on the contrary, cases with no oscillations of populations that are in strong coupling. 
The two concepts are therefore not equivalent, as none implies the other. This is illustrated 
for the SE case in Figs. [5^, [5]d and [6] on the one hand, where the system is in SC, and in 
Fig. [5t and [7] on the other hand, where it is in WC. In SS, there is no t dynamics in any 
case, so oscillations of populations are clearly unrelated to weak or strong coupling. In SE, 
the distinction is clearly seen in Fig. [5] where both the t and r dynamics are shown in a 
contour-plot in the case where the system is initially prepared as an exciton, (a) and (c), 
or as a polariton, (b). In the polariton case, the dynamics in t is simply decaying (because 
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T/g T/g r/g 

FIG. 5: (Color online) Time dynamics of the correlator ^{a^ {t)a{t + r))^^, cf. eqn (I27p . Only the 
pattern of oscillations is of interest here (lighter blues correspond to higher values). In all cases, 
both the t and r dynamics tend to zero. Figures (a) and (b) show the SE of an exciton and of an 
upper polariton, respectively, in a very strongly-coupled system (7^ = 0.2g and = O.l^i). Fig. (c) 
shows the SE of an exciton in weak-coupling (7^ = 5.9(7). The oscillations in r, rather than in t, 
are the mark of SC. 

of the lifetime), while it is clearly oscillating in r, were the proper manifestation of SC is 
to be found. The t decay is not exactly exponential because in the presence of dissipation, 
the polariton is not anymore an ideal eigenstate (the larger the dissipation, the more the 
departure). However this effect in SC is so small that it only consists in a small "wobbling" 
of the T contour lines. On the other hand, the exciton, (a), that is not an eigenstate, features 
oscillations both in the t dynamics (the one often but unduly regarded as the signature of 
SC), as well as the r dynamics. In stark contrast, the exciton in WC, (c), bounces with t. 
This, that might appear as an oscillation, is not, as it happens only once and is damped 
in the long-time values. This behaviour is shown quantitatively in Fig. [6] for SC and Fig. [7] 
for WC, where the population na{t) is displayed for the SE of an exciton (blue solid), a 
photon (purple dashed) and an upper polariton (brown dotted), respectively, along with the 
luminescence spectrum that they produce (detected in the cavity emission). Here it is better 
seen how, for instance, the polariton-decay is wobbling as a result of the dissipation, that 
perturbs its eigenstate-character and leaks some population to the lower polariton. More 
importantly, note how very different the spectra are, depending on whether the initial state 
is a photon or an exciton, despite the fact that the dynamics is similar in both cases (see the 
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inset in log-scale of their respective populations). The PL spectrum observed in the cavity 
emission is much better resolved when the system is initially in a photon state, than it is 
when the system is initially in an exciton state. The splitting is larger and the overlap of the 
peaks smaller in the former case. This will find an important counterpart in the SS case. In 
Fig. [71 the corresponding case of WC is shown for clarity, with a decay of populations and 
possible oscillations. 




FIG. 6: (Color online) Strong-Coupling spectra Sq^{lo) (left) and its corresponding mean number 
dynamics n^^{t) (right) for the SE of three different initial states: In blue solid, one exciton; in 
purple dashed, one photon and in brown dotted, one upper polariton. Parameters are 7a = 1.9g 
and 7fe = 0-lg. Inset of (b) is the same in log-scale. 

Figure [8] shows the r dynamics in the SS (when the t dynamics has converged and is 
steady), for five cases of interest to be discussed later (in Fig. [9]). A first look at the 
dynamics would seem to gather together a group of two curves that decay exponentially 
to good approximation (and remain positive as a result), and another group of three that 
assume a local minimum. The correct classification is the most counter-intuitive in this 
regard, as it puts together the dashed lines on the one hand and the solid on the other, 
i.e., scrambling them together. The mathematical reason for this classification is revealed 
in the inset, where the same dynamics is plotted on log-scale. The dashed (resp. solid) lines 
correspond to parameters where the system is in WC (resp. SC) according to our definition, 
i.e., to values of R that are imaginary on the one hand and real on the other. In log-scale, this 
corresponds respectively to a damping of the correlator, against oscillations with an infinite 
number of local minima. Note that the blue dashed line features one local minimum, which 
does not correspond to an oscillatory — or coherent-exchange — behaviour of the fields, but 
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FIG. 7: (Color online) Weak-Coupling spectra Sq^{uj) (upper) and its corresponding mean number 
dynamics n^{t) (lower) for the SE of an exciton (a) and a photon (b) as the initial condition. 
In all figures, solid blue corresponds to the decay of an exciton; dashed purple to the decay of 
a photon and dotted brown to the decay of an upper polariton. For comparison, we plotted in 
(a) and (b), with solid black lines, the very diff'erent bare emission [g = 0) of an exciton and 
photon respectively. Also in (a), in dashed black, that of an exciton decaying with the Purcell rate 
7p = 4(7^ /7a. Parameters are 7a = l.Qg and 7^ = O.lg. Inset of (c) is the same in log-scale. 

rather to a jolt in the damping. These considerations that may appear abstract at this level 
will later turn out to show up as the actual emergence of splitted (dressed) states or not in 
the emitted spectrum. 

We now return to the general (SE/SS) expression for the spectra, eqn fl30|) . that, at 
resonance in SC, simplifies to: 



So{uj) 



1 

2^ 



Rq 



{CI 



Kq 



(38) 
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FIG. 8: (Color online) Dynamics of limt_oo(aH*)a(i + T))/nf^, eqn ([27]) and ([Ha]), for the SS 
corresponding to the points (a)-(e) in Fig. [9j In inset, the same in log-scale. Solid lines (b, blue; c, 
purple and d, brown) of SC feature oscillations of the correlator, as the mark of SC. Dashed lines 
(a, green and e, blue) correspond to WC. Note that although the blue dashed line (e) appears to 
be similar to other SC lines, it does not oscillate in the log-scale, where it only features a single 
local minimum. In the same way, the brown line (d) that seems not to oscillate actually features 
an infinite set of local minima, as is revealed in the log scale. 

where we used the definition for the Rabi frequency at resonance, eqn ( 12^ . and 



In the weak coupling regime, with Rq pure imaginary {g < |r_|), the positions of the two 




(39b) 



(39a) 



peaks collapse onto the centre, ooa = = 0. Defining iR„ = Rq, with R„ = a/T^ 
real number, the general expression for the spectra rewrites as: 




1 9^{Do} 



(40) 



2tt R. 
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with the Lorentzian and dispersive contributions now given by: 



UJ 



We now address the specifics of the SE and SS cases. 



A. Case of Spontaneous Emission 

In the most general case of SE, the D^^ coefficient at resonance, -Dq^, is a complex 
number. If the initial condition further fulfils 3?n°j = 0, it becomes pure imaginary. Usually, 
the initial states considered are independent states of photons or excitons (not a quantum 
superposition), where indeed n^^ = [12I, Q. In these cases. 



which yields the following expression for the spectrum: 



nSE _ , 2V '^'^b i_b^ . 



Sri^) = -7 ' \ ' • (43) 



As noted earlier, the SE spectrum of exciton observed in the leaky modes is obtained 
from eqn ( H3ll by exchanging the indices a ^ b. We illustrate this with the two particular 
cases that follow. 

The typical detection geometry for the spontaneous emission of an atom in a cavity 
consists in having the atom in its excited state as the initial condition, and observing its 
direct emission spectrum. In this case the role of the cavity is merely to affect the dynamics of 
its relaxation, that is oscillatory with the light-field in the case of SC. This case corresponds 
to n° = 1 and n° = = in eqn (H3il with a ^ b. This gives [l^ : 

S^'^iu) = , ^ 4 +^ J _ (44) 

TT ^4 + ^2(7^ _ 2g^) + (g^ + 7a^)(^?2 + 2^)2 



In the semiconductor case, one would typically still have in mind the excited state of the 
exciton as the initial condition, but this time, this is the cavity emission that is probed. 
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The initial condition is therefore the same as before but without interchanging a and b in 
eqn fj43|) . which reads in this case: 

CSE/ \ _ 

The difference in the hneshape due to the initial quantum state is seen in Fig. [61 The 
visibility of the line-splitting is much reduced in the case of an exciton in SC which SE is 
detected through the cavity emission, than in the case of a photon. With a polariton as an 
initial state, only one line is produced. 

Again, by reason of symmetry, interchanging a ^ 6 in eqns (jH]) and (l45l) . correspond to 
the SE of the system prepared as a photon at the initial time and detected in, respectively, 
the cavity emission on the one hand (eqn HU a ^ b), and in the leaky mode emission on 
the other hand (eqn H5l) . In the latter case, the spectrum is invariant under the exchange 
a ^ b. Fig. ([6]) also hints to the changes brought by the detection channel (direct emission 
of the exciton or through the cavity mode). 

If n° = or 77,^ = (in which case n^^ = 0), the normalized spectra do not depend on the 
nonzero value n° or n°. That is, one cannot distinguish in the lineshape, the decay of one 
exciton from that of two, or more. In the more general case, when ra^^, 7^ 0, the peaks can be 
differently weighted. For instance, starting with an upper 

polariton 1 0o) = (|1, 0) + |0, l))/v^ 
(n° = n° = = 1/2) gives rise to a dominant upper-polariton peak (labelled 2 in the above 
equations, as seen in the brown dotted line in Fig. [6]). One can classify the possible lineshapes 
obtained for various initial states. For instance, as we have just mentioned, the normalized 
spectrum of \0,n) as an initial state, is the same whatever the nonzero n, which is not 
entirely unexpected from a linear model. From the previous statement, the same spectrum 
is also obtained for a coherent state or a thermal state of photons, or indeed any quantum 
state, as long as the exciton population remains zero. In the same way, the PL spectrum 
of the product of coherent states in the photon and exciton fields, \z) \z') with z = z' E C*, 
is the same as that of a polariton state |0, 1) + |1,0), although both are very different in 
character: a classical state on the one hand and a maximally entangled quantum state on 
the other. 
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B. Case of continuous, incoherent pumping 



In the SS, at resonance, Dq^ is pure imaginary: 

"\^(Pa + A)+W+' ^^^^ 

and the term that consists in the difference of Lorentzians in eqn ( l30ll disappears: ^J{C} = 0. 
As a result, the two peaks are equally weighted for any combination of parameters: 



(4' + CD + ^- (Al - A!) 

Kq 



(47) 



The only way to weight more one of the peaks than the other in the SS of an incoherent 
pumping, would be to pump directly the polariton (dressed) states, as is the case in higher- 
dimensional systems were polaritons states with nonzero momentum relax into the ground 



state 251] or in OD case when cross pumping is considered 2l|]. In our present model, 
however, such terms are excluded. The two peaks of the Rabi doublet, composed of a 
Lorentzian and a dispersive part, are both symmetric with respect to Ua = ujb = 0. Only 
if 53{-Dq^} = r_/(7, the spectrum of eqn (1471) consists exclusively of two Lorentzians. The 
parameters that correspond to this case are those fulfilling either = p^°p r^r^ or r_|_ = 0. 
The second case corresponds to the limiting case of diverging populations, where the SC 
becomes arbitrarily good. In the most general case, the dispersive part will contribute to 
the fine quantitative structure of the spectrum, bringing closer or further apart the maxima 
and thus altering the apparent magnitude of the Rabi splitting. In some extreme cases, 
as we shall discuss, it even contrives to blur the resolution of the two peaks and a single 
peak results, even though the modes split in energy. As for the weak coupling formula, it 



simplifies to: 



2ir V R, 

(48) 



2n V 

Both decompositions, eqns (H7j) and (l48l) . have been given to spell-out the structure of the 
spectra in both regimes. The unified expression that covers them both reads explicitly: 
CSS. ^ ^ J 8^^n + 2P.(4^^ + rg) 

° nnf IQiu^ - 4u;2(8^2 _ r2 _ + (4^2 + r^r^)2 • i^^^ 

It is the counterpart for SS of eqn (143|) . for SE. The case of excitonic emission can also be 
obtained, as for SE, by interchanging the indices a ^ b. 
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C. Discussion 



Although the spectra in the semiconductor case that are probed at neghgible electronic 
pumping [Ph ^ 1) with no cavity pumping at all {Pa = 0), are in principle described by the 
same expression as that of the SE case used in the atomic model, in practise, however, both of 
these conditions can be easily violated. The renormalization of 7^ with P^ brings significant 
corrections well into the linear regime. For instance, for parameters of point (c) in Fig. [9] 
with Pa = 0, the rate Pb that is needed to bring a 100% correction to 7;, yields, according 
to eqns fll4ap and (I14bp . averages much below unity, namely, n^^ ^ 0.026 and ~ 0.121. 
By the time Uf, reaches unity, with Ua still one fourth smaller, the correction on the effective 
decay rate has became 400%. As this is Ua which is proportional to the signal detected 
in the laboratory, the electronic pumping must be kept very small so that corrections to 
the effective linewidth can be safely neglected. As the more interesting nonlinear regime is 
probed, the renormalized Fs become very different from the bare 7s. 

Second, even in the vanishing electronic pumping limit, it must be held true that Pa 
is zero. Even if only an electronic pumping is supplied externally by the experiment, the 
pumping rates of the model are the effective excitation rates of the cavity and exciton field 
inside the cavity, and it is clear that photons get injected in the cavity in structures that 
consists of numerous spectator dots surrounding the one in SC (cf. Fig. [2]). Although most 
of these dots are in WC with the cavity, they affect the dynamics of the SC QD by pouring 
cavity photons in the system. In the steady state, following our previous discussion, this 
corresponds to changing the effective quantum state for the emission of the strongly-coupled 
QD. As we shall see in more detail in what follows, this bears huge consequences on the 
appearance of the emitted doublet, especially on its visibility. 

To fully appreciate the importance and deep consequences of these two provisions made by 
the SS case on its SE counterpart, we devote the rest of this Section to a vivid representation 
in the space of pumping and decay rates parameters. Now that it has been made clear what 
is the relationship of the SE case with respect to the SS one, we shall focus on the latter 
that is the adequate, general formalism to describe SC of QDs in microcavities. 

In presence of a continuous, incoherent pumping, the criterion for strong-coupling — from 
the requirement of energy splitting and oscillations in the r dynamics that we have discussed 
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above — gets upgraded from its usual expression [12l |: 

^7>l7-l, (50) 

to the more general condition: 

g>\T.\. (51) 

The quantitative and qualitative implications and their extent are shown in Fig. [9], where 
we have fixed the parameters 7^ = O.lg and Pa/g = 0.5g, and outlined the various regions 
of interest as Ph and 7^ are varied. This choice of representation allows us to investigate 
the changes that can be easily imprinted experimentally in the system: by tuning on a 
system of varying quality with regard to SC, as specified by its quality factor that tunes 7a. 

The red lines, that enclose the filled regions, delimit a frontier above which the pump is so 
high that populations diverge (there is no steady state). This is expressed mathematically 
by the condition that the common denominator of the mean values in eqns ( !T4|) vanishes. It 
is equivalent to setting the total decay rate of the system — that should remain positive in 
order for the correlators to decay in time — to zero. In the SC regime, such a rate is given 
by r+ and therefore the condition for a SS reads: 

Pb = la + 7b-Pa, (i.e., r+ = 0). (52) 

In the WC regime, the system decays with an effective Purcell rate and the condition trans- 
forms into 

= ^ ^ p . (53) 

7a + 7b - Pa 

The main separation inside that region where a SS exists, is that between SC (in shades 
of blue, inside the triangle) and WC (in shades of red, on its right elbow). The blue solid 
line that marks this boundary, specified by 

Pa = ^g-la + lb + Pa, (i.e. , = | | ) , (54) 
separates the weak and strong coupling regions. The dashed vertical black line, specified by 

7a = % + 7fe, (i-e., ^ = I7-I) , (55) 

corresponds to the standard criterion of SC (without incoherent pumping). 

The light-blue region, labelled 1 on Fig. [9l corresponds to SC as it is generally under- 
stood. It has a clear splitting of the lines in the luminescence spectrum. The dark-blue 
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region, labelled 2, corresponds to SC, according to the requisite that Rq be real, but with a 
broadening of the lines so large that in the luminescence spectrum, eqn (j47j), only one peak 
is resolved. This region is delimited by the brown line, which is the solution of the equation 
(PS{uj)/duj'^\i^=Q = 0, i.e., with no concavity at the origin. From this condition follows the 
implicit equation: 

(3r+ - r_)/ + (r_ - r+f + g\Do\ig' - r^^ + sr^) = o (56) 

that yields two solutions, only one of which is physical (the other is placed on the red 
line r+ = 0, where the system diverges). Note that this line extends into the WC region, 
as we shall discuss promptly. The distinction between line-splitting as it results from the 
emergence of new dressed states in the SC, and the observation of two peaks in the spectrum, 
is seen clearly on Fig. [10], where the two are superimposed and seen to differ greatly even at 
a qualitative level for most of the range of parameters, coinciding only in a narrow region. 
The doublet as observed in the luminescence spectrum collapse much before SC is lost. Any 
estimation of system parameters, such as the coupling strength, from a naive association 
of the PL spectrum with that of two-coupled oscillators, will most likely be off by a large 
amount. 

The last region of SC, labelled 3, is that specified by Ag + ■yi, < ■ya < 4:g + % + Pa — Pb, 
i.e., that which satisfies eqn ( l5Ti) but violates eqn ( 150|) . thereby being in SC according to 
the more general definition that takes into account the effect of the incoherent pumping, 
but that, according to the conventional criterion, is in WC. For this reason, we refer to this 
region as of pump-aided strong coupling. This is a region of strong qualitative modification 
of the system, that should be in WC according to the system parameters, but that restores 
SC thanks to the cavity photons forced into the system. 

We now consider the other side of the blue line, that displays the counterpart behaviours 
in the WC. Region I is that of WC in its most natural expression. Region II, in light, is the 
extension into WC of the characteristic of featuring two maxima in the emission spectrum. 
In this case, this does not correspond to a line-splitting in the sense of SC where each peak 
is assigned to a renormalized (dressed) state, but rather to a resonance of the Fano type 
that is carving a hole in the single line of the weakly-coupled system. In this region, one 
needs to be cautious not to read SC after the presence of two peaks at resonance. Finally, 
region III is the counterpart of 3, in the sense that this region is in SC according to the 
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conventional criterion, eqn (150|) . for the system parameters, when in reahty the too-high 
electronic pumping has bleached the strong-coupling. 

In inset, we reproduce the diagram to position the five points (a)-(e) in the various 
regions discussed, for which the luminescence spectra are displayed and decomposed into 
their Lorentzian and dispersive contributions, eqns (1391) . Case (c), at the lower- left angle, 
corresponds to SC without any pathology nor surprise: the doublet in the luminescence 
spectrum — although displaced in position as shown in Fig. [10] — is a faithful representation 
of the underlying Rabi-splitting. Increasing pumping brings the system into region 2 where, 
albeit still in SC, it does not feature a doublet anymore. The reason why, is clear on the 
corresponding decomposition of the spectrum. Fig. MJo), with a broadening of the dressed 
states (in purple) too large as compared to their splitting. Further increasing the pump 
brings it out of the SC region to reach point (a), where the two Lorentzians have collapsed. 
This degeneracy of the mode emission means that the coupling only affects perturbatively 
each mode. As a result, the dispersive correction has vanished, and the spectrum now 
decomposes into two new Lorentzians centred at zero, with opposite signs. 

Back to point (c), now keeping the pump constant and increasing 7^, we reach point (d). 
It is still in SC, although the cavity dissipation is very large (more than four times the 
coupling strength) for the small value of 7;, considered. Its spectrum of emission shows, 
however, a clear line-splitting that is made neatly visible thanks to the cavity (residual) 
pumping Pa- Increasing further the dissipation eventually brings the system into WC, but 
in region II where, again due to Pa 7^ 0, the spectrum remains a doublet. In Fig. (e), one 
can see, however, that there is no Rabi splitting, and that the two peaks arise as a result of 
a subtraction of the two Lorentzians centred at zero [see the WC spectrum decomposition 
in eqn (140114 lal) ]. There is probably no need to display a spectrum from region I, as in this 
case it does not show any qualitative difference as compared to that of (a). 

The set of system parameters and the estimated effective pumping rates allow to re- 
construct Fig. [9] through eqn (15214561) . that contains all the physics of the system. In the 
following, we shall look at variations of this representation to clarify or illustrate those 
aspects that have been amply discussed before. 

Fig. [11] shows the same diagram as that of Fig. [9], only with Pa now set to zero, i.e., 
corresponding to the case of a very clean sample with no spurious QDs other than the SC 
coupled one, that experiences only an electronic pumping. Observe how, as a consequence. 
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region 3 of SC and II of WC have disappeared. The former was indeed the result of the 
residual cavity photons helping SC. The "pathology" in WC of featuring two peaks at reso- 
nance in WC has also disappeared, but most importantly, see how region 2 has considerably 
increased inside the "triangle" of strong-coupling, meaning that the parameters required 
so that the line-splitting can still be resolved in the luminescence spectrum now put much 
higher demands on the quality of the structure. This difficulty, especially in the region 
where Pb g follows from the "effective quantum state in the steady state" , that we have 
already discussed. The presence of a cavity pumping, even if it is so small that no field- 
intensity effects are accounted for, can favour SC by making it visible, indeed by merely 
providing a photon-like character to the quantum state. This is the manifestation in a SS 
of the same influence that was observed in the SE: the luminescence spectrum of a photon 
as an initial state of the coupled system, is more visible than that of an exciton, keeping 
all parameters otherwise the same (see Fig. [6]). Another useful picture to highlight this last 
point, is that where the various regions are plotted in terms of the pumping rates. Pa and Pb 
[see Fig. [121 for the case 7a = 3.8g of the line (a)-(c) in Fig. ([9])]. The angle of a given point 
with the horizontal, linked to = Pb/Pa, deflnes the exciton-like or photon-like character 
of the SS established in the system, and thus determines the visibility of the double-peak 
structure of SC. This is, at low pumpings, independent of the magnitude \/ Pa + Pbi ^ 
the brown line deflned by eqn fl56l) is approximately linear in this region. This shows the 
importance of a careful determination of the quantum state that is established in the SS by 
the interplay of the pumping and decay rates, through eqn ( |T2l) . The magnitude, on the 
other hand, affects the splitting i?o, and the linewidths r+. In order to have a noticeable 
renormalization, the pumps must be comparable to the decays. On the one hand, the Rabi 
frequency can be affected in different ways by the pumpings, depending on the parameters. 
If Ta = Tb, there is, in general, no effect of decoherence on the splitting, showing that in this 
case there is a perfect symmetric coupling of the modes into the new eigenstates (although 
the broadening can be large and spoil the resolution of the Rabi splitting anyway). If they 
are different, for example in the common situation that 7a — 76 > Pa — Pb, the Rabi increases 
with increasing Pa — Pb- On the other hand, the linewidth F.,. = {'~fa + lb — Pa — Pb)/4: presents 
clear bosonic characteristics: it increases with the decays but narrows with pumping 
The intensity of the pumps also affects the total intensity of the spectra, that is proportional 
to n^^ through 7a and the integration time of the apparatus. Here, however, we have focused 
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on the normalised spectra (i.e., the hneshape). 

V. SUMMARY AND CONCLUSIONS 

In conclusion, we have brought under a unified formalism the zero-dimensional light- 
matter interaction, both in the Weak (WC) and Strong (SC) Coupling, for the two cases 
of Spontaneous Emission (SE) of an initial state, and emission under a Steady State (SS) 
maintained by an incoherent, continuous pumping. While the SE case for some particular 
initial conditions (excited state of the atom) and configurations (resonance, direct emission) 
has been a pillar of SC in cQED of atoms in cavities, the extensions that we provided 
here to include a continuous and incoherent pumping, are suitable to describe the recent 
field of cQED in semiconductor heterostructures. Together, they merge into an elegant and 
complementary theoretical edifice. 

The main results of this papers are to be found in eqns (I30l)-(l33l) that provide the an- 
alytical expression for the cavity emission spectra of a system whose specificities — such as 
whether it corresponds to SE or the SS established by an incoherent continuous pumping — 
are provided by a parameter which, in the first (SE) case, is given by eqn flMl) . and in the 
second (SS), by eqn (!36|) . These formulas, that allow for an arbitrary detuning between the 
bare modes, reduce to more self-contained expressions at resonance, namely eqn fl43|] for SE 
and eqn (H9l) for SS. The resonance case allows an unambiguous definition of SC, depending 
on whether the complex Rabi frequency, eqn (l28l) . is pure imaginary (WC) or real (SC). This 
corresponds, in turn, to a damping or to sustained oscillations of the time-autocorrelation 
of the fields. This is completely independent of the dynamics of the populations. SC is 
characterised by the emergence of new eigenstates, with different energies, whereas in WC, 
the energies remain degenerate. There is no, however, one-to-one mapping of this splitting 
of the energies with the lines observed in the luminescence spectrum. All cases can arise: 
one or two peaks can be observed at resonance both in WC and SC. In the SC case, one 
peak only is observed when the energy splitting is too small as compared to the broadening 
of the lines, whereas in the WC, two peaks are seen as a result of a resonance carving a 
hole in a single line, giving the illusion of a doublet (and indeed of an anticrossing when 
detuning is varied). For that reason, an understanding of the general picture is required to 
be able to position a particular experiment in the space of parameters, as was done in Fig. [9] 
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and Figs. [TTHT2I rather than to rely on a quahtative effect of anticrossing at resonance. 
Figure [To] shows how loosely related are the observed line-splitting in the luminescence spec- 
trum (solid red) and the actual energy splitting of the polariton modes (dressed states, in 
dotted blue). The various situations that may arise are illustrated and discussed in Fig. O 
The respective effects of the angle and the distance to origin in the P^, Pf, parameter-space, 
is shown in Fig. [121 the angle accounts for the effective quantum state, that imparts on the 
visibility of the splitting in the spectrum, while the magnitude accounts for bosonic effects 
like line- narrowing and field-intensity renormalization of the Rabi splitting. 

This work addresses the case of bosonic excitons, that corresponds to the case of vanishing 
excitations, even for genuine two-level excitons. As such, it also contains a lot of the physics 
of ground state quantum wells excitons in planar cavities, although in this case, another type 
of pumping — a so-called cross-terms pumping — is more adequate, as particles are injected 
directly into the ground state by scattering of polaritons. On the contrary, in our present 
scheme, the excitation is in terms of the bare modes, through phonon-assisted scattering of 
electron-hole pairs into the QD for the electronic pumping, Pb, or via Purcell-emission of 
weakly-coupled spectator dots into the cavity mode for the cavity pumping. Pa- A natural 
extension of this work is to consider fermionic excitons, that do not admit more than one 
particle. In this case, the equations of motions for the correlators are not closed, and only 
semi-analytical results are available. The structure of the spectra — that in the bosonic case 
decompose as a Lorentzian and a dispersive line for each peak — becomes that of an infinite 
series of lines, tightly grouped together, to give rise to multiplet structures with no more 
than four peaks, albeit in a wide variety of different configurations. This fermionic case 
corresponds to mapping the Jaynes-Cummings ladder to an exact luminescence spectrum, 
in much the same way that we have been doing with the Rabi doublet in this text. This 
case that is of crucial importance for the study of nonlinearity of genuine (two-levels) QDs 
in semiconductor microcavities, is the subject of part II of this text. 
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APPENDIX A: TIME-RESOLVED DYNAMICS IN THE SE CASE 



From eqn (fT2!) . the mean photon population for the SE of a system whose initial conditions 
are given by eqn f[T^ . is: 



n^(t) = e-2^+*<J cos {V2zit) + cosh (v^Za^) 
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The expression for ni,{t) follows from a ^ b. The crossed mean value that reflects the 
coherent coupling reads: 
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where we have defined the complex parameters: 



z^,2=\l\R\'±ig'-i' 
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In case of SC, at resonance, {g > |7~|), both zi^2 are real, due to the fact that > |7-A| 
for whatever the detuning, and they can be written in a more transparent way: 
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The mean value na{t) is plotted in Fig. [6](b) and[7](c) for SC and WC, respectively, for the 
parameters in the caption. 

It might be of interest to note that eqns (lAlHASp . but also eqn fl27l) in the SE case, and 
therefore all the results that follow from them, are reproduced by introducing decay as an 
imaginary part to the energies in the Heisenberg picture, i.e., upgrading Ua^h to uja,b — ila,b/'^ 
and solving directly in a full Hamiltonian picture the operator equations of motion: dtc{t) = 
i[H,c{t)] with c = a,b. This simple expedient resurges plainly for instance in eqn fl28l) . 
that is an exact and important result derived with the full dissipative formalism, but which 
follows directly from eqn ([6]) with complex energies. This recourse to complex energies is 
however not valid in the Schrodinger picture and/or with the pumping terms. 



APPENDIX B: DERIVATION OF THE LUMINESCENCE SPECTRUM IN THE 
SS FROM THE GENERAL EXPRESSION 

Expression (fT6|) is general, but its final form in term of a time-integrated Fourier trans- 
forms of G^^\ eqn f|T9|l . is not rigorous for the SS case since both the numerator and the 
denominator are infinite quantities. Their ratio, however, produces a finite quantity, which 
bears all the characteristic of the optical emission spectrum in a SS, and — may be more 
convincingly — recovers the well-established and rigorously derived Wiener-Khint chine for- 
mula, eqn (120|) . Note that this formula is, strictly speaking, an arcane mathematical result in 
the theory of stochastic processes. There, S(uj) is a mesure of the strength of the fluctuations 



of the Fourier component at frequency u 19|. It has no strict connection with a physical 



signal, as both infinite negative and positive times are required for its demonstration, which 
violates causality among other things. For a rigorous and extended discussion of a physical 



optical spectrum, see the excellent discussion by Eberly and Wodkiezicw [26|. As the step 
from the general eqn f|T6|) to the SS, eqn fl20l) . can be made straightforwardly, we provide it 
below: it merely consists in cancelling the diverging quantities by obtaining the final result 
as a limit of the time-integrated spectrum. We believe that the starting point being so clear 
mathematically, and this presentation so satisfying in its simultaneous coverage of both the 
SE and SS cases, it is worth considering for a pedagogical purpose. 

The population and the autocorrelator can be decomposed as a transient (TR) and steady- 
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state (SS) values: 

{a^a){t) = (a'^a)™(t) + lim {a^a){t) , (Bla) 

i— »oo 

, r) = {a^t)a{t + r))™ + lim {a^{t)a{t + r)) , (Bib) 

t—*oo 

where limt^oo('3^a)(t) = n^^, cf. eqn f ll4al) . We rewrite eqn f|T9l) as the time integration of 
the Fourier transform until time T, that is left to increase without bounds: 

' {a^a) {t)dt Jt=o Jt=o 

Substituting eqns (IBip in this expression, we can keep track of the terms that cancel (one 
can check, from the explicit results of the text, the convergence of the quantities = 
/Q°'(ata)™(t)rft and = J^{a^{t)a{t + r))^^-- dt, for all r > 0): 
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Since the norm of the Fourier transform of is also bounded (that, again, can be checked 
from the explicit result), the limit in T yields: 



1 

^ lim 3? / {a\t)a{t + T))e"^^dT , (B4) 
K Jo 

which is the formula of the text, eqn fl20|) . 
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u)lg u/g u/g 

FIG. 9: (Color online) Phase space of the SS Strong/Weak Coupling as a function of Pb/g and 
^a/g for the parameters 76 = Q.lg and Pa/g = 0.5g. The red lines, eqns ([52H53]) . delimit the region 
where there is a steady state. The blue line, eqn (j54p . separates the strong (in shades of blue) from 
the weak (shades of red) coupling regions. The dotted black line, 7a = 4g' + 7^, separates SC and 
WC regions in the absence of pumping. The brown line, eqn (j56p . separates the regions where one 
(dark blue) or two (light blue) peaks can be resolved in the luminescence spectra. This delineates 
three areas in the SC region: (1) two peaks are resolved in the spectra, (2) the two peaks cannot be 
resolved and effectively merge into one, albeit in SC, and (3) SC is achieved thanks to the pump Pa 
(with one or two peaks visible depending of the overlap with the light or dark area) despite the 
large dissipation that predicts WC according to eqn (jSOp . In the same way we can distinguish three 
regions in weak coupling: (I) standard WC, (II) SC with a two peaked spectrum and (III) WC 
due to pumping P^. The surrounding figures (a) to (e) show spectra from these regions and their 
decomposition into, Lorentzian and dispersive parts in SC, or Lorentzians centred at zero in the 
WC. Parameters correspond to the points in the inset: (a) 7a = 3. 85 and P}, = 5, (b) 7a = i.Sg and 
Pb = 0.55, (c) 7a = 3.89 and n = O.lg, (d) 7a = 4.2^ and A = O.I5, (e) 7a = 4.8^ and P^ = O.lff. 
Observe how, in SC, two eigenstates have emerged in purple dotted lines, even in the cases — like 
in (b) — where they are not seen in the total spectrum. In the same way, in WC, all the emission 
emanates from the origin, although a two-peak structure can arise as a result of a resonance, also 
centred at the origin. 
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FIG. 10: (Color online) Rabi slitting at resonance (dotted blue) ±iio) eqn (f29|) . and the observed 
position of the peaks in the PL spectra (solid red) as a function of Pb/g- Parameters are those of the 
line of points (a), (b) and (c) of Fig. O 7a = 3.8(7, lb = 0.1(/, Pa = 0.5g, and the corresponding Pb. 
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FIG. 11: (Color online) The same phase space of SC/WC as in Fig. [9]as a function of Pb/g and 7a/(7, 
only with Pa/g now set to zero (no cavity pumping). The triangle of strong coupling is displaced, 
and the regions of one peak spectra (regions 2 of SC and I of WC) are enlarged, following eqn ()56p . 
As a result regions 3 and II have disappeared. 
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FIG. 12: (Color online) Phase space of SC/WC as a function of the pumps Ph/g and Pa/ 9 for fixed 
decay parameters 7a = 3. 85 and 75 = O.lg. As in Fig. [9l the red lines mark the physical regions 
and the blue one the SC (blue shades)/WC (red shades) transition, with the same regions 1 and 2 
of SC and III of WC, also with the points (a), (b) and (c), of Fig. [9j In inset, zoom of the low-pump 
region, showing the importance of both the angle, arctan (1/a), and the magnitude, \J Pa ~^ of 
a given point. 
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